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ABSTRACT 



(N 



There is an ongoing debate in cosmology about the value of the length scale at which 
"homogeneity" in the matter distribution is reached or even if such a scale exists. In 
the wake of this debate, we intend in this letter to clarify the meaning of the statement 
■ transition to homogeneity and of the concept of correlation length. We show that there 

are two scales each associated to the fractal and to the homogeneous regimes of the 
. matter distribution, respectively. The distinction between both scales has deep conse- 

quences, for there can be a regime which despite having small fluctuations around the 
OO 

\ average density exhibits large clusters of galaxies. 

Subject headings: methods: statistical — galaxies: clusters: general — large-scale struc- 
+3 , ture of universe 



1. Introduction 

The analysis of the distribution of large-scale structure in the Universe is done by applying 
statistical tools. The most basic descriptor is, aside from the average mass density, p = (/j(r)), the 
reduced two-point correlation function ( Peebles 1993[ ), 



(p(r)p(O)) (5p(r)5p(0)) 
£( r ) = J2 1 = J2 ' C 1 ) 

where we assume statistical homogeneity and isotropy. Higher order correlation functions of the 
probability distribution are also of importance, but their determination becomes less reliable as the 
order increases. However, for the purpose of illustrating our point, it will suffice to concentrate 
on £(r). The study of galaxy catalogues leads to fitting £(r) by a power law on the length scales 
probed by these catalogues, 

««=(?) 7 . P> 
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but the values of tq and 7, and even whether the latter is scale-dependent, are still being de- 
bated (Coleman & Pietronero 1992; Sylos Labini, Montuori Pietronero 1998| ; Davis 1997; Guzzo 
1997; |Scaramella et al. 1998| ; Martinez et al. 1998| ). At any rate, it seems that scale-invariance is 
approximately true on a non-negligible range of scales. 

The distance ro, defined by S,( r o) = 1 (or by some other qualitatively equivalent condition, 
e.g. |Martmez et al. 1998| ), has b een called "correlation length" in the literature. However, it is 



well known in the theory of critical phenomena ( Stanley 1971 ; Binney et al. 1992 ; Ma 1994 ) that 
power-law correlated fluctuations let themselves be felt over the entire range of observation. There, 
the physical meaning of a correlation length corresponds to the length scale where power-law decay 
crosses over to exponential decay, and has nothing to do with the condition £ ~ 1. 

The distance ro has certainly the meaning of separating a regime of large fluctuations, 5p 3> p, 
from a regime of small fluctuations, bp <C p. The regime of small fluctuations has been identified 
with a homogeneous distribution and, because of this, ro has been called the scale of transition 
to homogeneity. However, as we will see, the sense in which bp <C p leads to homogeneity is not 
equivalent to absence of structure. In fact, fluctuations correlated by a power law, albeit small, may 
give rise to collective, macroscopic phenomena, such as for instance critical opalescence in fluids 
(Stanley 1971; Binney et al. 1992] ) . Absence of structure is not determined just by the condition 
that correlations vanish at large distances, but rather by how fast they do so. This is measured by 
the correlation length Ao, separating slow from fast decay. Fast decay is generally identified with 
exponential decay, but it could also refer to another type of sufficiently rapid decay. 

The point, therefore, is to establish a clear distinction between the strength of the correlations 
(a measure of how large the fluctuations around the mean can be) , and the range of the correlations 
(a measure of how far correlations are felt). These two characteristics of the correlation function, 
though related, are not equivalent to each other. This subtle difference has been usually overlooked 
in the cosmological literature or, at least, not clearly stated. 



2. Relevant concepts 

In order to exemplify the main concepts, we consider a random distribution of particles whose 
average density is p. In the next section we will discuss the implications for cosmology. The second 
moment of the density field is given by 

G{r) = p 2 [l + i{r)]. (3) 

Correlations are weak if £ < 1, and strong if £ > 1. Thus the scale ro defined by £(ro) = 1, 
which separates these two regimes, is related to the strength of the correlations. In the parlance 
of structure formation, ro is the scale separating at a given time the linear from the nonlinear 
dynamical regimes (ISahni fe Coles 1995| ), and we may call it the scale of nonlinearity. A more 



intuitive picture can be achieved if one examines the variance of the smoothed density field, pl, the 
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number of particles contained in a volume of size L 3 per unit volume. This variance is computed as 
an integral of £(r). We consider the cosmologically interesting case where the correlation function 
has a simple scale- invariant form like (|2|), with given r$ and 7 < 3 (this latter requirement assures 
convergence of the integral). Then we can write: 

where we have defined a length scale a = /cP 1 / 3 (a 3 is the volume per particle), and B is a positive 
numerical constant of order one. The first term is the shot-noise contribution, and the second is due 
to correlations. We conclude that if the volume containing the sample is large enough [L 3> f*o,a), 
the relative fluctuations of pi around its mean are very small, and a good estimate of p can be 
extracted from the sample. What happens for intermediate values of L depends on the relative 
values of a and Tq. We consider two physically interesting cases: (i) a < ro, as may be the case 
with the galaxy distribution, so that for scales a <C L <^ ro fluctuations are large and the correlation 
term dominates; and (ii) a ~ ro, as in a fluid in thermal equilibrium, so that for scales L <^ a, ro, 
fluctuations are large but dominated by shot noise. In either case no reliable estimate of p can be 
extracted. We also see that a Gaussian approximation to the one-point probability distribution 
of pl is certainly excluded in these cases, because pl > by definition and hence the probability 
distribution becomes prominently skewed if 8pL ^ P- 

On the other hand, the range of the correlations is related to how fast they decay at infinity, 
which is measured by the spatial moments of the two-point correlation function: 

M n = ( drr rf (r) oc lim V k V k P(k), (5) 

where we have introduced the Fourier transform of the two-point correlation function, the power 
spectrum P(k) = J dr e* kr £(r). Notice that A4 n vanishes for odd n due to isotropy of the correla- 
tions. Starting from the interpretation of £(r) as the density contrast of particles at a distance r 
from any one particle, it is straightforward to give a physical meaning to the moments A4 n : they 
are the multipoles characterizing the morphology of the typical cluster generated by correlations. 
Correlations are short-ranged if all the moments are finite, which implies that they decay faster 
than any power of the distance and that P(k) is analytic at k = 0. One can then define a typical 
size of the cluster, Aq, the correlation length, as 



\ 2 



1 


TrM 2 


= lim — 


VlP(k) 


2 


M 


k^o 2 


P(k) 
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in agreement with the standard definition in condensed matter physics (see e.g. Ma 1994] ) ; if P(0) 
or V 2 i P(0) vanishes, then the definition of Ao must be generalized in terms of a different quotient 
between higher-order, non-vanishing moments. Therefore, if one observes a realization of the field 
ql with a smoothing length L S> Ao, then clusters are unobservable and the distribution is indistin- 
guishable from the case that there are no correlations at all. Conversely, if £(r) decays as a power 
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of r, so that it is long-ranged, then the moments M n diverge from some n onwards (depending on 
the exponent of decay), and P(k) is not analytic at k = 0. Hence, no matter how large L is, the 
existence of clusters will always be detectable, so that one associates an infinite correlation length 
to a power-law decay. In conclusion, it is Ao the scale that marks the transition to homogeneity in 
the sense of absence of structure. Notice that in principle Ao and ro are different and in fact they 
can even be orders-of-magnitude different, as a power-law decay of £(r) exemplifies. 

It is clear now that there are in general at least three different length scales (a, ro, Ao), thus 
defining four regimes, which we now discuss. Notice that without further knowledge of the physics 
of the system, one cannot exclude that any two of the scales are of the same order, which is not a 
rare case in physics, so that some of the regimes below can be difficult to observe. In what follows 
we will assume the case most appropriate in a cosmological context, that a < r^ < Xq. We then 
have: 



Homogeneous regime, for scales r ^> Ao- On these scales, fluctuations are small (£(r) <C 1) 
and the system structureless. Typically, correlations are exponentially damped and pi has 
a prominently peaked distribution, so that the mean density p can be reliably estimated. 
An example is a liquid in thermal equilibrium far from a phase transition, for which £(r) ~ 
(ro/r) exp(— r/Ao), where a, ro and Ao are all of the order of the molecular diameter ( [Landau] 
fc Lifshitz 1980| ; |Goodstein 19851 ). 

Critical regime, for scales ro r Ao- Fluctuations are still small, so that the mean 
density p can be estimated with confidence too but now because fluctuations are correlated 
and therefore behave coherently on these scales, they give rise to spatially extended struc- 
tures. Moreover, if the correlations follow a power law on these scales, fluctuations about the 
background have a self-similar, fractal structure. An example is a fluid at the critical point, 
for which £(r) is like in equation (|2|) with 7 < 3. Again a and ro are of the order of the 
molecular diameter, but now Ao — > 00. Notice that small fluctuations do not imply uninter- 
esting behavior: since they are macroscopic in spatial extent, one observes phenomena such 
as critical opalescence ( [Stanley 1971 ; Binney et al. 19921) . The reader may feel uneasy with 



the apparent contradiction between the claim that fluctuations are small in a critical fluid and 
the standard claim that they are in fact large. The point is to identify with which quantity 
one is comparing the fluctuations. When compared with p, as we do, they are small, which 
is a prerequisite to render a thermodynamic description valid. When compared with typical 
fluctuations far from the critical point, they are extremely large, and in fact the ratio (size of 
critical fluctuations) / (size of non-critical fluctuations) usually diverges in the thermodynamic 
limit, as is seen from equation (||). 

Fractal regime, for scales a <C r <C ro. Now fluctuations exhibit structure and at the same 
time are very large (£ 3> 1), so that from equation (]3|) it follows that G(r) ~ p 2 £,{r), and 
p cannot be estimated. For power-law correlations, this means that it is the density itself 
(and not only its fluctuations around the mean) that appears to have fractal structure (or 
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multifractal, depending on how higher-order correlations scale). This manifests itself in that 
most realizations of the stochastic field p(r) do not look like a continuum field at all on these 
scales, but are very irregular and full of voids of every size separated by very high density 
condensations. This regime is not encountered in a critical fluid because ro is of the order of 
the molecular diameter, but may be relevant in cosmology, as we will see. 

• Shot-noise regime, for scales r <C a. The presence of correlations is masked by shot noise: 
fluctuations are very large, but simply because the discrete nature of the underlying point- 
process is dominant. 

It may be helpful at this point to show pictorially a distribution in the critical regime. Figure 1 
is a snapshot of a critical two-dimensional lattice gas,[] with Ao — > oo, but a and ro are both of the 
order of a pixel (=atom), barely discernible to the naked eye. In a critical fluid the average density 
is well defined, but one can clearly see in the figure structures (i.e., inhomogeneities) on scales 
much larger than ro. Although one may observe what seem like voids of arbitrary size, a closer 
examination unveils that the voids contain clusters of atoms. A fractal however must contain voids 
of every size, owing to scale invariance, and the voids need to be totally empty (density exactly 
zero in them), so the fluctuations of the density over a volume of size L are large at any scale. It is 
now clear that what has scale invariance, or fractal character, in a critical fluid is the fluctuations 
over the average density, rather than the density field itself. For comparison, Figure 2 shows the 
same model at very high temperature, with Ao of the order of a pixel (comparable to a and ro), and 
where it exhibits Poissonian fluctuations. Visual inspection shows immediately that fluctuations 
are less prominent than in the critical regime, and that there is no structure at all. 



We now apply the above to cosmology. We have two different scales, ro and Ao, related to the 
intensity and to the range of the correlations, respectively. This does not preclude the existence of 
other intermediate length scales, which do not however affect the physical meaning of ro and Ao, nor 
the discussion to follow. Incidentally, one of these scales, which has also been termed correlation 



where <f> is the peculiar gravitational potential due to the fluctuations in the density field, 5p. This 
length can be understood as the characteristic scale of spatial variation of the field (p. But this 



3. Implications for cosmology 



length, is defined as (Kofman et al. 1993; Sathyaprakash et al. 1995[ ) 




(7) 



It actually is isomorphic to the two-dimensional Ising model. The three-dimensional Ising model as a description 
of the galaxy distribution has been introduced in Hochberg & Perez- Mercader 1996; Perez-Mercader et al. 1996. 
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is not the correlation length either, because a correlation between the values of the field at two 
separate points does not imply in general that these values must be similar. 

As discussed in the previous section, the length scale Ao is given by the behavior of the power 
spectrum on the largest scales. Since one cannot observationally gather information about scales 
larger than the horizon, Ao can be estimated only by assuming that the known behavior of P(k) on 
the largest observable scales holds beyond the horizon. The power spectrum on the largest scales 



can be extracted from the analysis of the microwave background, yielding P(k) ~ k (|Bennett et 



ah 1996 ; 5 moot 1999 ) at the epoch of decoupling. If this non-analytic behavior were extrapolated 
to k — > 0, it would imply Ao = oo. The physical interpretation of this conclusion is that the 
correlation length was certainly larger than the horizon at the epoch of decoupling. Hence, since 
causal processes can create correlations up to the scale of the horizon only, we must be observing the 
remaining, frozen imprints of a pre-inflationary epoch which have not been erased yet. Moreover, 
given that the evolution of the power spectrum by gravitational instability does not alter its shape 
in the linear regime ( Peebles 1993| ), one must conclude that Aq is still beyond the horizon in the 



present epoch. Thus, it may seem that the homogeneous regime is beyond reach and that the largest 
observable scales in the Universe belong to the critical regime. As a matter of fact, notice that the 
well-known pictures of the microwave background by COBE ( Bennett et al. 1996|) resemble figure 1 
much more than figure 2. 

Let us consider now the scale of nonlinearity, ro- The smallness of density fluctuations as 
inferred from observations of the microwave background ( {(5p L ) 2 )/p 2 L ~ 1(T 7 for L ~ 10 3 h' 1 
Mpc) yields an upper bound on ro. Smaller scales can be probed by means of the galaxy catalogues. 
The conclusions related to the matter distribution depend on the bias parameter, but, in agreement 
with theoretical models, this parameter can be assumed to be of order unity. This implies in turn 
that ro can be directly extracted from the galaxy-galaxy correlations. The standard result therefore 
is that ro ~ 5 h^ 1 Mpc ( pavis fc Peebles 19831 ; Peebles 1993j ), so that currently available galaxy 



catalogues, of size L S> ro, are probing the critical regime. Hence, the galaxy distribution should 
look very much like a fluid at the critical point: small but spatially extended fluctuations. As 
already stressed, this does not mean that the catalogues are probing the transition to homogeneity, 
since ro is not the correlation length. 

There has recently arisen, however, a controversy on the validity of ro ~ 5 h^ 1 Mpc, due to 
suggestions that ro is in fact larger than the maximum effective size of the current largest galaxy 
catalogues, and that a simple power-law form like (0) with a fixed exponent holds at least up to this 
maximum scale ( poleman Pietronero 1992j ; |Sylos Labini, Montuori fc Pietronero 19"98| ) . If so, 
then the catalogues are probing the fractal regime. As shown in the previous section, this implies 
that the mean galaxy number density cannot be reliably estimated from these catalogues, and hence 
neither can ro. This means that they do not allow to discriminate between a distribution that has a 
nonvanishing p and looks fractal up to scales of the order of ro , and a distribution whose p vanishes 
and thus is fractal up to arbitrarily large scales. Not surprisingly, the idea of a universe with a 
vanishing average density as the size of the sample volume increases is very old, dating back to at 
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least 1908 ( |Charlier 1908|) . This construction was based on a hierarchy of clustering without limit, 



that is, clusters of galaxies which form clusters of clusters, which in turn also form clusters and so 
on, ad infinitum. It is illustrative to read now the old review by G. de Vaucouleurs advocating for 
a hierarchical cosmology (de Vaucouleurs 1970). 



As we have seen before, it appears that the correlation length Ao is much larger than the 
scale of nonlinearity, tq. This could lead to interesting consequences. As already remarked, and 
illustrated by the discussion of figures 1 and 2, there can exist structure in the critical regime, i.e. 
on intermediate scales ro -C r <C Xq. Hence, this could be the reason why clustering is observed, as 
galaxy clusters and superclusters, even on scales much larger than ro, which is a "paradox" in the 
interpretation of galaxy catalogues. Similarly, the increase of ro with bias threshold should not be 



interpreted in terms of a larger correlation length for the corresponding structure (Gabrielli, Sylos 



Labini, & Durrer 1999) 



Another consequence of a large Aq is related to the effect of gravitational lensing on the 



microwave background (|Blanchard & Schneider 1987 


; Kashlinsky 1988; 


Cole & Efstathiou 198$ 


Martinez-Gonzalez, Sanz, & Silk 1990; 


Fukushige, Makino, & Ebisuzaki 1994 


; Martinez-Gonzalez, 


Sanz, & Cayon 1997 


). The geodesic equation for light propagation is formally identical to the 



geometric-optics equation for light-rays in a medium of non-uniform refractive index, which allows 
one to identify the gravitational potential with an effective refractive index. Therefore, the fact that 
the correlation length of the potential diverges means that gravitational lensing of the microwave 
background could exhibit an effect akin to critical opalescence (a sort of "cosmological opalescence" , 
Dommguez & Gaite 1999p . 



In conclusion, we have provided a clear explanation and a precise mathematical formulation 
of the difference between the intensity of fluctuations and their spatial extent in the context of 
cosmological structures. These two properties can be characterized by two different length scales: 
the scale of nonlinearity, ro, and the correlation length, Ao- We have discussed the behavior of the 
correlations in the different regimes which these two lengths define and the connection with cosmo- 
logical observations. It should be clear now that the statement "transition to homogeneity" as used 
in the cosmological literature really means "transition to the regime of small-amplitude fluctua- 
tions", which does not imply absence of structure. In fact, we have also argued that this remaining 
"critical structure" could have non-trivial consequences for the interpretation of observations. For 
instance, no real transition to homogeneity is likely to be observed in the galaxy distribution. 
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Fig. 2: Non-critical fluid. 



